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Semi-metallic  Be5C2  monolayer  global  minimum 
with  quasi-planar  pentacoordinate  carbons  and 
negative  Poisson's  ratio 


Yu  Wang\  Feng  i\\  Yafei  Li^  &  Zhongfang  Chen^ 


Designing  new  materials  with  novel  topological  properties  and  reduced  dimensionality  is 
always  desirable  for  material  innovation.  Here  we  report  the  design  of  a  two-dimensional 
material,  namely  86502  monolayer  on  the  basis  of  density  functional  theory  computations.  In 
86502  monolayer,  each  carbon  atom  binds  with  five  beryllium  atoms  in  almost  the  same 
plane,  forming  a  quasi-planar  pentacoordinate  carbon  moiety.  86502  monolayer  appears  to 
have  good  stability  as  revealed  by  its  moderate  cohesive  energy,  positive  phonon  modes  and 
high  melting  point.  It  is  the  lowest-energy  structure  with  the  8e502  stoichiometry  in 
two-dimensional  space  and  therefore  holds  some  promise  to  be  realized  experimentally. 
8e502  monolayer  is  a  gapless  semiconductor  with  a  Dirac-like  point  in  the  band  structure  and 
also  has  an  unusual  negative  Poisson's  ratio.  If  synthesized,  8e5C2  monolayer  may  find 
applications  in  electronics  and  mechanics. 
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Carbon  in  known  molecules  and  materials  typically  has 
tetrahedral  tetracoordination  (for  example,  diamond), 
planar  tricoordination  (for  example,  graphite)  or  linear 
dicoordination  (for  example,  ethyne)  arrangements.  Especially, 
the  tetrahedral  preference  of  tetracoordinate  carbon  compounds, 
which  was  deduced  by  van’t  Hoff^  and  Lebel^  more  than  a 
century  ago,  is  one  of  the  major  foundations  of  organic 
chemistry.  In  1968,  Monkhorst^  first  discussed  an  exceptional 
bonding  pattern  of  carbon,  namely  planar  tetracoordinate 
carbon  (ptC),  by  proposing  a  planar  methane  (D4h),  which  is 
actually  not  a  minimum  energy  structure.  Later  in  1970,  by 
insightful  analysis  of  ptC  bonding  in  the  hypothetical  planar 
methane,  Hoffman  et  al.^  suggested  that  ptCs  can  be  stabilized 
electronically  by  replacing  H  atoms  in  planar  methane  by 
G  donors  (to  facilitate  electron  transfer  to  electron -deficient  a 
bonds)  or  7i  acceptors  (to  delocalize  the  unfavourable  lone  pair  of 
ptC).  Along  this  line,  in  1976  Collins  et  al^  theoretically  devised 
the  first  ptC-containing  molecule  1,  1-dilithiocyclopropane.  One 
year  later  Cotton  and  Millar^  synthesized  the  first  ptC-containing 
molecule,  although  the  unique  ptC  configuration  was  not 
recognized  at  that  time.  Ever  since,  ptC  chemistry  has  been  a 
subject  of  extensive  studies^”^.  Numerous  ptC  species  have  been 
designed  theoretically^ and  some  global  minimum  structures, 
such  as  CAI4+  (ref.  16),  CAI4+  (ref.  17)  and  CAl3Si“  (ref.  18) 
have  been  observed  experimentally.  More  excitingly,  besides  ptC, 
many  molecules  containing  planar  pentacoordinate  carbon 
(ppC)^^“^^  and  hexacoordinate  carbon  (phC)^^“^^  have  been 
designed  computationally.  The  rule-breaking  chemical  bonding 
in  these  planar  hypercoordinate  carbons  can  lead  to  completely 
new  molecules  and  materials,  which  are  of  fundamental 
importance  to  chemistry  and  materials  science. 

The  unique  topology  and  dimensionality  may  lead  to 
exceptional  properties  of  materials.  It  is  not  a  surprise  to  witness 
the  growing  interest  in  designing  ptC-containing  solids 
and  nanostructures  in  recent  years^^“^®.  Especially,  stimulated 
by  the  extensive  studies  of  graphene^  and  inorganic  layered 
materials^^’^^,  many  two-dimensional  (2D)  materials  with 
rule-breaking  chemical  bonding  have  been  designed 
computationally^^”"^^.  Eor  example,  on  the  basis  of  ptC 
molecule  CB4  (ref.  43),  Wu  et  al?^  designed  the  first 
ptC- containing  2D  material,  namely  B2C  graphene,  which 
was  later  confirmed  to  be  only  a  local  minimum^^.  Li  et 
and  Dai  et  al?^  predicted  that  AI2C  monolayer  in  its 
lowest-energy  configuration  have  all  C  atoms  being  ideal  ptC. 
More  interestingly,  Li  et  demonstrated  recently  that  one 
C  atom  can  bind  to  six  beryllium  (Be)  atoms  in  an  almost  planar 
manner,  yielding  a  phC-featuring  2D  material  with  a  Be2C 
stoichiometry.  With  intriguing  structural  and  electronic 
properties,  these  2D  materials  are  expected  to  have  important 
applications  in  some  specific  fields. 

In  contrast  to  ptC  and  phC,  no  attempt  exists  for  extending 
ppC  into  solids  until  now.  Motivated  by  the  successful  design  of 
ptC-  and  phC-containing  materials,  we  are  quite  curious  whether 
it  is  possible  to  extend  ppC  into  solids.  Addressing  this  issue 
would  deliver  not  only  unique  structures  but  also  some  fantastic 
properties,  which  is  of  both  theoretical  and  practical  importance. 

In  this  work,  by  means  of  density  functional  theory  (DET) 
computations,  we  first  confirm  that  the  ppC  molecule  Be9C2  “  is 
a  local  minimum  on  the  potential  energy  surface.  Based  on  the 
structural  characters  of  Be9C2  “ ,  we  design  a  ppC-containing  2D 
material,  namely  Be5C2  monolayer,  in  which  each  C  atom  binds 
to  five  Be  atoms  to  form  a  quasi-planar  pentacoordinate  moiety. 
Our  Be5C2  monolayer  has  rather  good  thermal  and  kinetic 
stabilities  and  is  energetically  the  most  favourable  isomer  in  2D 
space.  Dramatically,  the  electronic  band  structure  of  Be5C2 
monolayer  has  a  Dirac-like  point  at  the  Eermi  level,  endowing  an 


intriguing  semimetallic  feature.  More  interestingly,  Be5C2 
monolayer  possesses  a  negative  Poisson’s  ratio,  which  is  rather 
rare  in  nanostructures.  These  fascinating  properties  make  Be5C2 
monolayer  a  promising  candidate  for  future  applications  in 
electronics  and  mechanics. 


Results 

Be9C2~  as  the  inspiring  ppC  species  for  2D  monolayer.  Our 

design  of  periodic  2D  Be5C2  monolayer  is  initially  inspired  by  our 
finding  of  the  ppC  molecule,  namely  Be9C2  ~ ,  which  has  singlet 
ground  state  and  D2h  symmetry  (Eig.  la).  Computed  at  B3LYP 
level  of  theory  with  6-311  +  ^^(d,  p)  basis  set,  Be9C2  ~  is  a  local 
minimum  with  the  lowest  vibrational  frequency  of  158.5  cm 
In  this  molecule,  each  C  atom  binds  to  five  Be  atoms  to  form  a 
ppC  moiety  of  BesC,  which  has  been  shown  to  be  a  local 
minima^^.  Structurally,  Be9C2“  can  be  viewed  as  two  BesC 
moieties  fused  by  sharing  one  Be  atom.  The  C-Be  bond  lengths 
are  in  the  range  of  1.70-1.76  A,  whereas  the  Be-Be  bond  lengths 
are  in  the  range  of  2.02-2.39  A. 

According  to  the  natural  population  analysis,  each  ppC  of 
Be9C2“  possesses  a  2.18  |e|  negative  charge  (  —  0.40  |e|  according 
to  the  Hirshfeld  charge  population  analysis)  and  the  natural 
electron  configuration  is  2s^-^^2pj^-^^2p];-^^2p^  indicating  that 
ppCs  in  Be9C2  are  essentially  stabilized  through  Be  cr-donation 
and  delocalization  of  carbon  electrons.  The  computed  Wiberg 
bond  index  for  C-Be  bonds  are  0.52  and  0.61,  respectively, 
resulting  in  a  total  Wiberg  bond  index  of  2.81  for  each 
ppC.  Moreover,  we  also  scrutinized  the  molecular  orbitals  of 
Be9C2  ~  to  get  more  information  on  its  stabilization  mechanism. 
As  shown  in  Pig.  lb,  the  highly  delocalized  n  (for  example, 
HOMO-7  and  HOMO-8)  and  g  (for  example,  HOMO-3, 
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Figure  1  |  BegCf  “  molecule,  (a)  Optimized  geometric  structure  of 
BegCf  “  molecule.  The  black  and  green  balls  represent  C  and  Be  atoms, 
respectively.  The  red  and  black  numbers  are  the  Wiberg  bond  index  (WBI) 
and  lengths  (in  angstroms)  of  the  representative  chemical  bonds, 
respectively,  (b)  The  canonical  molecular  orbitals  of  BegCf”. 
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HOMO -4  and  HOMO-5)  orbitals  can  help  maintain  the  planar 
configuration. 


Geometric  properties  of  Be5C2  monolayer.  The  fusing  of  two 
BesC  moieties  to  Be9C2~  reminds  us  of  the  roadmap  of  fusing 
benzene  rings  to  polycyclic  aromatic  hydrocarbons  (for  example, 
naphthalene  and  anthracene)  and  then  to  2D  infinite  graphene. 
Inspired  by  the  ppC-containing  Be9C2  “ ,  we  designed  a  new  2D 
inorganic  material,  namely  Be5C2  monolayer  by  generalizing  the 
structural  characters  of  Be9C2  “ .  As  shown  in  Fig.  2a,  1  unit  cell 
of  Be5C2  monolayer  consists  of  8  C  atoms  and  20  Be  atoms  with 
the  optimized  lattice  parameters  being  a  =  8.92  A  and  ^  =  9.21  A, 
respectively.  Similar  to  Be9C2  “ ,  in  Be5C2  monolayer  each  C  atom 
binds  to  five  Be  atoms  to  form  a  ppC  moiety  of  BesC  and  two 
neighbouring  BesC  moieties  share  one  Be  atom  (BeJ  to  form  a 
Be9C2  moiety  in  the  a  (Be9C2-I)  or  b  (Be9C2-II)  direction. 
Especially,  one  Be9C2-I  moiety  is  fused  with  four  neighbouring 
Be9C2-II  moieties  by  sharing  the  peripheral  Be  atoms  (Be2)  and 
vice  versa,  leading  to  the  formation  of  a  2D  network  with  four 
Be9C2  moieties  in  one  unit  cell.  The  optimized  coordinates  of 
Be5C2  monolayer  are  presented  in  Supplementary  Table  1.  It  is 
noteworthy  that  the  Be5C2  monolayer  is  remarkably  buckled 
rather  than  an  exactly  planar  structure  (Fig.  2b).  The  buckling, 
measured  by  the  vertical  distance  between  the  bottommost  Be 
atoms  and  the  uppermost  Be  atom,  is  as  high  as  2.14  A.  Even  so, 
the  ppCs  in  Be5C2  monolayer  still  have  a  good  planarity 
(Supplementary  Fig.  1).  According  to  our  computations,  the  total 
degrees  of  five  Be-C-Be  angles  for  ppCs  in  Be9C2-I  and  Be9C2-II 
moieties  are  371.44°  and  364.84°,  respectively,  which  are  a  little 
higher  than  the  ideal  360°.  Interestingly,  when  Be9C2~  is 
neutralized  by  four  protons  (H  +  ),  the  obtained  Be9C2H4 
molecule  is  also  severely  buckled  (Supplementary  Fig.  2).  In 
Be5C2  monolayer,  the  length  of  C-Bei  bonds  (1.66  A,  1.68  A) 
is  much  shorter  than  that  of  C-Be2  bonds  (1.73  A,  1.74  A). 
Moreover,  our  computations  revealed  that  Be5C2  monolayer  has  a 
non-magnetic  ground  state,  indicating  that  there  are  no  unpaired 
electrons  in  Be5C2  monolayer. 

We  then  computed  the  deformation  electronic  density  of  Be5C2 
monolayer  to  elucidate  its  bonding  nature.  The  deformation 
electronic  density  is  defined  as  the  total  electronic  density 
excluding  those  of  isolated  atoms.  As  clearly  shown  in  Fig.  2c, 

a  c 


Figure  2  |  Geometric  and  electronic  structures  of  60502  monolayer. 

(a)  Top  and  (b)  side  views  of  geometric  structure  of  Be5C2  monolayer.  The 
blue  dashed  lines  denote  a  unit  cell;  a  and  b  represent  the  lattice  vectors,  1 
and  2  denote  the  different  Be  atoms,  I  and  II  represent  the  different  Be9C2 
moieties,  (c)  Deformation  charge  density  of  2D  Be5C2  monolayer.  Red  and 
yellow  refer  to  electron  accumulation  and  depletion  regions,  respectively. 


some  electrons  are  extracted  from  the  2s  orbitals  of  Be  atoms  and 
well  delocalized  over  C-Be  bonds,  indicating  that  C  atoms  form 
multicentre  covalent  bonds  with  neighbouring  Be  atoms,  which  is 
crucial  for  stabilizing  the  ppC  moieties.  The  similar  stabilizing 
mechanism  has  been  found  in  ptC-  and  phC-containing  2D 
materials^^”"^^’"^^.  According  to  the  Hirshfeld  charge  popular 
analysis,  C,  Bei  and  Be2  atoms  in  Be5C2  monolayer  possess  a 
—  0.32,  +0.16  and  +0.12  |e|  charge,  respectively.  The  buckling 
of  Be5C2  monolayer  stretches  the  Be-Be  distances  and  probably 
helps  reduce  the  otherwise  even  stronger  repulsive  interactions 
between  Be  atoms. 

We  also  used  the  recently  developed  Solid  State  Adaptive 
Natural  Density  Partitioning  (SSAdNDP)  method^^  to  better 
understand  the  unique  chemical  bonding  of  Be5C2  monolayer. 
According  to  our  results  (Supplementary  Fig.  3),  there  is  no 
classical  two-centre-two-electron  (2c-2e)  C-Be  bond  in  Be5C2 
monolayer.  For  one  unit  cell  of  Be5C2  monolayer,  the  SSAdNDP 
search  revealed  twenty-four  3c-2e  Be-C-Be  cr-bonds  (responsible 
for  the  bonding  within  the  BesC  units),  four  4c-2e  cr-bonds  on 
four  Be  squares  and  eight  6c-2e  7i-bonds  over  eight  BesC  units, 
accounting  for  72  electrons  per  unit  cell.  This  bonding  pattern  is 
consistent  with  the  symmetry  of  Be5C2  monolayer.  Especially,  the 
existence  of  delocalized  a-  and  71 -bonds  could  essentially  stabilize 
the  ppCs  in  Be5C2  monolayer. 

Stability  of  Be5C2  monolayer.  Although  Be5C2  monolayer  has 
rather  intriguing  structural  properties,  we  are  unclear  whether  it  is 
a  stable  structure.  To  assess  the  stability,  we  first  computed  the 
cohesive  energy  of  Be5C2  monolayer,  which  is  defined  as: 
Fcoh  =  (^£c  +  ^£Be-£Be5C2)/(^  +  ^)>  in  which  Fo  ^Be  and 
^Be5C2  are  the  total  energies  of  a  single  C  atom,  a  single  Be  atom 
and  Be2C  monolayer,  respectively;  n  and  m  are  the  number  of  C 
and  Be  atoms  in  the  supercell,  respectively.  According  to  our 
computations,  Be5C2  monolayer  has  a  cohesive  energy  of  4.58  eV 
per  atom.  As  a  reference,  the  cohesive  energies  of  the  experi¬ 
mentally  realized  2D  silicene^^’^^  and  phosphorene"^^’^®  are  3.71 
and  3.61  eV  per  atom,  respectively.  As  silicene  and  phosphorene 
are  composed  of  covalent  bonds,  the  even  higher  cohesive  energy 
can  ensure  that  Be5C2  monolayer  is  a  strongly  connected  network. 

The  stability  of  Be5C2  monolayer  can  be  further  confirmed  by 
its  phonon  dispersion  curves.  As  shown  in  Fig.  3a,  no  appreciable 


Figure  3  |  Kinetic  stability  of  Be5C2  monolayer,  (a)  Phonon  spectrum  and 
(b)  phonon  DOS  of  80502  monolayer. 
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imaginary  phonon  mode  is  present,  suggesting  the  good  kinetic 
stability  of  Be5C2  monolayer.  Remarkably,  the  highest  frequency 
of  Be5C2  monolayer  reaches  up  to  1,120  cm  which  is 
higher  than  those  of  M0S2  monolayer  (473  cm  silicene 

(580  cm  and  our  recently  designed  phC-containing  Be2C 

monolayer  (1,020  cm  “  indicating  robust  C-Be  bonds  in 
Be5C2  monolayer.  The  analysis  of  partial  phonon  density  of  states 
(DOS;  Fig.  3b)  revealed  that  the  highest  frequency  of  Be5C2 
monolayer  is  mainly  contributed  by  C-Bei  bonds. 

Moreover,  to  examine  the  thermal  stability  of  Be5C2 
monolayer,  we  performed  first-principles  molecular  dynamic 
(FPMD)  simulations  using  a  2x2  supercell.  Our  three 
simulations  at  temperature  of  1,000,  1,500  and  2,000  K  show 
that  Be5C2  monolayer  can  maintain  its  structural  integrity 
throughout  a  10-ps  FPMD  simulation  up  to  1,500  K,  but  is 
seriously  disrupted  at  2,000  K,  suggesting  that  Be5C2  monolayer 
has  a  melting  point  between  1,500  and  2,000  K  (Supplementary 
Fig.  4).  We  also  performed  optimizations  for  the  distorted 
structures  from  MD  simulations.  After  full  atomic  relaxation, 
those  structures  from  MD  simulations  at  1,000  and  1,500  K  can 
recover  the  initial  configuration  of  Be5C2  monolayer.  The  above 
results  demonstrate  that  Be5C2  monolayer  has  a  remarkable 
thermal  stability. 

The  moderate  cohesive  energy,  all  positive  phonon  modes  and 
good  thermal  stability  can  ensure  that  Be5C2  monolayer  is  at  least 
a  local  minimum  structure  on  the  potential  energy  surface. 
However,  is  the  ppC-containing  Be5C2  monolayer  a  global 
minimum?  It  is  noteworthy  that  the  global  minimum  structure 
is  more  likely  to  be  achieved  than  the  local  minimum  structures 
experimentally.  For  example,  many  isomers  of  graphene,  such  as 
Haeckelite  graphene^^,  T-graphene^^  and  penta-graphene^^  have 
been  designed  computationally,  but  none  of  them  has  been 
realized  experimentally.  Therefore,  we  performed  a  global  search 
for  the  lowest-energy  structure  of  Be5C2  monolayer  in  the  2D 
space  using  first-principles-based  particle-swarm  optimization 
(PSO)  method  as  implemented  in  CALYPSO  code.  After  a 
comprehensive  search,  we  obtained  three  stable  isomers  of  2D 
Be5C2,  which  are  labelled  as  Be2C5-I,  Be2C5-II  and  Be2C5-III, 
respectively.  As  shown  in  Fig.  4a,  Be2C5-I  is  actually  the  above 
discussed  ppC-featuring  Be5C2  monolayer.  Interestingly,  in 
Be2C5-II  (Fig.  4b)  and  Be2C5-III  (Fig.  4c),  each  C  atom  binds  to 
five  Be  atoms  to  form  a  ppC  moiety  of  BesC  and  the  Be9C2 
moieties  also  can  be  found  in  these  two  isomers.  Be2C5-I  is  50  and 
101  meV  per  atom  lower  in  energy  than  Be2C5-II  and  Be2C5-III, 
respectively,  indicating  that  Be2C5-I  is  the  global  minimum 


a  b  c 


Figure  4  |  Low-energy  isomers.  Geometric  structures  of  80502-!  (a), 
80502-!!  (b)  and  80502-!!!  (c).  The  b!ue  dashed  !ines  denote  unit  ce!!s. 


structure  in  the  2D  space.  Therefore,  the  ppC-featuring  Be5C2 
monolayer  holds  great  potential  to  be  realized  experimentally. 

Electronic  properties  of  Be5C2  monolayer.  With  such  interesting 
structural  characteristics,  does  Be5C2  monolayer  also  have 
intriguing  properties?  To  address  this  issue,  we  computed  the 
band  structure  of  the  lowest-energy  Be5C2  monolayer.  As  shown 
in  Fig.  5a,  Be5C2  monolayer  is  gapless  or  semi-metallic  with  the 
conduction  band  minimum  (CBM)  and  valence  band  maximum 
(VBM)  meeting  at  the  Fermi  level,  which  is  quite  similar  to  that  of 
graphene.  However,  the  meeting  point  of  Be5C2  monolayer  is 
located  on  the  path  from  G  (0,  0,  0)  point  to  Y  (0,  0.5,  0)  point 
rather  than  on  a  high-symmetry  point  as  for  graphene.  Especially, 
the  conduction  and  valence  bands  around  the  Fermi  level  exhibit 
a  linear  dispersion  relation,  suggesting  that  the  meeting  point  of 
Be5C2  monolayer  is  also  Dirac-like.  Considering  that  the  PBE 
functional  tend  to  underestimate  the  band  gap,  we  recomputed 
the  band  structure  of  Be5C2  monolayer  using  the  hybrid 
HSE06  functional^^  and  found  that  the  dispersion  of  the  valence 
and  conduction  bands  at  the  Fermi  level  is  similar  to  that 
predicted  by  PBE  and  no  appreciable  band  gap  can  be  identified 
(Supplementary  Fig.  5).  Thus,  the  gapless  property  of  Be5C2 
monolayer  is  solid. 

To  obtain  deeper  insight  into  the  electronic  properties  of  Be5C2 
monolayer,  we  analysed  its  DOS.  As  shown  in  Fig.  5b,  the  DOS  is 
zero  at  the  Fermi  level,  which  further  supports  the  presence  of  the 
Dirac  cone.  The  partial  DOS  analysis  shows  that  the  VBM  and 
CBM  are  contributed  by  both  Be-2p  and  C-2p  states,  and  the 
contribution  from  Be-2p  states  is  much  more  than  that  from  C-2p 
states.  Moreover,  we  plotted  the  partial  charge  densities  of  the 
VBM  and  CBM.  As  shown  in  Fig.  5c,d,  both  VBM  and  CBM  are 


Figure  5  |  Electronic  properties  of  60502  monolayer,  (a)  8and  structure 
and  (b)  DOS  of  80502  mono!ayer.  The  Fermi  !eve!  is  assigned  at  OeV. 
(c,d)  The  isosurfaces  of  partia!  charge  densities  for  the  (c)  V8M  and 
(d)  C8I\/1  of  80502  monoiayer.  The  isovaiue  is  0.015 e  A 
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mainly  originated  from  the  delocalized  orbitals  of  Be  atoms  and 
partially  from  the  multicentre  bonding  between  C  and  Be  atoms. 

Mechanical  properties  of  Be5C2  monolayer.  Besides  the 
electronic  properties,  we  also  investigated  the  mechanical 
properties  of  Be5C2  monolayer  by  examining  its  elastic  constants. 
As  a  validation,  the  computed  elastic  constants  of  graphene  are 
Cl  1  =  C22  =  342.93  GPa  and  C12  =  C21  =  62.23  GPa  respectively, 
which  achieve  good  agreements  with  experimental  measure¬ 
ments^^  and  previous  computations^^.  For  Be5C2  monolayer,  its 
elastic  constants  were  computed  to  be  Cn  =  32.90  GPa, 
C22=  130.89  GPa,  Ci2  =  C2i=  -5.32GPa  and  Gee  =  48.32  GPa, 
which  are  in  agreement  with  the  mechanical  stability  criteria  for  a 
tetragonal  2D  sheet  (C11C22  —  Ci2>0,  C66>0)^^.  The  in-plane 
Young’s  modules  along  a  (Y^)  and  b  (Y^)  directions,  which  can  be 
deduced  from  the  elastic  constants  by  Y^  =  (C11C22  —  ^12^21)1^22 
and  Y^  =  (C11C22  —  Ci2C2i)/Cii,  were  computed  to  be  32.68 
and  130.03  Nm“  L  respectively.  As  Y^  is  not  equal  to  Y^,  Be5C2 
monolayer  is  mechanically  anisotropic.  Moreover,  computed  at 
the  same  level  of  theory,  the  in-plane  Young’s  modules  of 
Be5C2  monolayer  are  less  than  those  of  graphene 
{Ya  =  Y^  =  331.63  N  m  “  but  higher  than  those  of 
phosphorene  (Y^  =  25.50 Nm“  ^  and  Y^  =  91.61  Nm“  ^), 
suggesting  that  Be5C2  monolayer  has  good  mechanical  properties. 

Remarkably,  we  noted  that  Be5C2  monolayer  has  a  negative 
C12,  which  results  in  negative  Poisson’s  ratios  of  —  0.041 
(C12/C22)  and  —0.16  (C12/C11)  for  a  and  h  directions, 
respectively.  It  is  noteworthy  that  the  Poisson’s  ratio  is  defined 
as  the  negative  ratio  of  transverse  to  axial  strain.  The  negative 
Poisson’s  ratio  indicates  that  Be5C2  monolayer  can  be  compressed 
or  stretched  in  both  two  directions  at  the  same  time.  For  a 
validation,  we  applied  a  uniaxial  strain  of  5%  in  a  and  h  directions 
of  Be5C2  monolayer,  respectively.  Just  as  expected,  the 
equilibrium  lattice  parameters  of  h  and  a  directions  are  elongated 
by  ~0.2  and  ~0.8%,  respectively,  confirming  that  Be5C2 
monolayer  indeed  has  negative  Poisson’s  ratios. 

Discussion 

The  unusual  negative  Poisson’s  ratio  may  endow  Be5C2 
monolayer  with  enhanced  toughness  and  shear  resistance, 
as  well  as  enhanced  sound  and  vibration  adsorption. 
Correspondingly,  Be5C2  monolayer  could  find  some  important 
applications  in  the  fields  of  mechanics,  tissue  engineering  and 
national  security.  It  is  worth  noting  that  the  negative  Poisson’s 
ratio  is  rather  peculiar,  as  in  nature  nearly  all  materials  have  a 
positive  Poisson’s  ratio,  except  some  so-called  auxetic  materials^^. 
Recently,  Jiang  et  demonstrated  theoretically  that  single¬ 
layer  black  phosphorus  has  a  negative  Poisson’s  ratio  due  to  the 
unique  puckered  configuration.  However,  the  negative  Poisson’s 
ratio  of  phosphorene  was  observed  in  the  out-of-plane  direction, 
which  is  different  from  that  of  Be5C2  monolayer.  Remarkably,  the 
Poisson’s  ratio  of  Be5C2  monolayer  in  the  b  direction  (  —  0.16)  is 
much  higher  than  that  of  phosphorene  ( —  0.027)^®,  rendering 
Be5C2  monolayer  a  more  promising  candidate  for  specific 
application  in  mechanical  devices.  No  doubt,  the  negative 
Poisson’s  ratio  of  Be5C2  monolayer  should  be  originated  from 
its  intriguing  structural  properties,  especially  the  uniquely 
oriented  chemical  bonds.  Our  results  could  provide  some 
guidelines  for  designing  materials  with  a  negative  Poisson ’s  ratio. 

With  so  many  fascinating  properties,  it  is  desirable  to 
synthesize  Be5C2  monolayer  in  the  laboratory.  Considering  that 
there  is  no  layered  structure  of  Be5C2  in  nature,  a  promising 
approach  is  to  grow  Be5C2  monolayer  on  the  surface  of  metal  or 
metal  oxide  via  chemical  vapour  deposition  with  accurately 
controlled  Be/C  ratio,  just  similar  to  the  growth  of  silicene^^’^^.  It 


is  noteworthy  that  Be  has  toxic  properties  and  the  chemical 
vapour  deposition  synthesis  usually  requires  high  temperature; 
hence,  special  caution  should  be  given  during  the  experimental 
realization. 

To  summarize,  inspired  by  the  bonding  characters  of  a 
ppC-containing  molecule,  Be9C2~,  we  designed  a  ppC- 
containing  2D  inorganic  material  with  a  Be5C2  stoichiometry. 
Our  DFT  computations  demonstrated  that  ppCs  in  Be5C2 
monolayer  are  essentially  stabilized  by  the  charge  transfer  from 
Be  ligands.  The  moderate  cohesive  energy,  absence  of  imaginary 
modes  in  phonon  spectrum  and  high  melting  point  evaluated 
from  FPMD  simulations  indicated  that  Be5C2  monolayer  is 
experimentally  viable.  Especially,  a  global  minimum  search 
revealed  that  Be5C2  monolayer  is  the  lowest-energy  structure 
for  the  Be5C2  stoichiometry  in  2D  space,  which  endows  Be5C2 
monolayer  great  possibility  to  be  realized  experimentally. 
Our  computations  demonstrated  that  Be5C2  monolayer  is 
semi-metallic  with  a  zero  band  gap  in  the  electronic  band 
structure.  More  interestingly,  Be5C2  monolayer  has  rather 
intriguing  mechanical  properties  featured  with  a  negative 
Poisson’s  ratio.  Therefore,  Be5C2  monolayer  is  expected  to  have 
wide  applications  in  electronics  and  mechanics.  We  hope  our 
theoretical  studies  will  promote  the  experimental  realization  of 
this  novel  material  and  attract  more  attentions  on  investigating 
nanomaterials  with  novel  chemical  bonding. 

Methods 

DFT  computations.  For  the  BegC^”  molecule,  geometry  optimizations,  frequency 
analyses  and  electronic  structure  computations  were  performed  at  the  B3LYP^^’^^ 
level  of  theory  with  the  6-311  -f  G^{d,  p)  basis  set  as  implemented  in  Gaussian  03 
package^^.  For  2D  Be5C2  monolayer,  DFT  computations  were  performed  using  the 
plane-wave  technique  implemented  in  Vienna  ab  initio  simulation  package^^.  The 
ion-electron  interaction  was  described  using  the  projector- augmented  plane  wave 
approach^^’^^.  The  generalized  gradient  approximation  expressed  by  PBE 
functional^^  and  a  500-eV  cutoff  for  the  plane-wave  basis  set  were  adopted  in  all 
the  computations.  The  convergence  threshold  was  set  as  10  “^eV  in  energy  and 
10  “  ^  eV  A  “  ^  in  force.  We  set  the  x  and  y  directions  parallel  and  the  z  direction 
perpendicular  to  the  layer  plane,  and  adopted  a  supercell  length  of  15  A  in  the  z 
direction.  The  Brillouin  zones  was  sampled  with  an  8  x  6  x  1  T  centred  k  points 
grid.  The  phonon  spectrum  was  computed  using  finite  displacement  method  as 
implemented  in  CASTEP  code^^.  The  elastic  constants  were  also  computed  using 
the  CASTEP  code.  The  chemical  bonding  analysis  of  Be5C2  monolayer  was  done 
using  the  SSAdNDP  method^^,  which  can  well  interpret  the  chemical  bonding  in 
terms  of  classical  lone  pairs,  two-centre  bonds,  as  well  as  multi-centre  delocalized 
bonds  in  bulk  solids,  surfaces  and  nanostructures. 

Molecular  dynamics  simulations.  The  thermal  stability  of  Be5C2  monolayer  was 
evaluated  by  means  of  first-principles  molecular  dynamics  simulations.  The 
ground-state  structure  of  Be2C  monolayer  was  annealed  at  different  temperatures. 
At  each  temperature,  MD  simulation  in  NVT  ensemble  lasts  for  10  ps  with  a  time 
step  of  1.0  fs.  The  temperature  was  controlled  by  using  the  Nose-Hoover  method^^. 

Global  minimum  structure  searches.  The  PSO  method  within  the  evolutionary 
scheme  as  implemented  in  the  CALYPSO  code^®  was  employed  to  find  the  low- 
energy  structures  of  2D  Be5C2  monolayer.  In  our  PSO  simulation,  the  number  of 
generation  was  maintained  at  30.  Unit  cells  containing  total  atoms  of  7,  14  and  28 
were  considered.  The  structure  relaxations  during  the  PSO  simulation  were 
performed  using  Vienna  ab  initio  simulation  package  at  PBE  level  of  theory. 

Data  availability.  All  relevant  data  are  available  from  the  authors. 
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